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Abstract 

The quadrupolar ordering in CeBg is explained in terms of the electrostatic 
interaction of quadrupolar moments arranged into a simple cubic lattice. The 
representation of magnetic and quadrupolar moments by means of quasispins 
of two kinds is employed. A linear increase of the quadrupolar transition tem- 
perature Tq{H) with applied magnetic field and its further re-entrance are 
described using a generalized spherical model which is well adjusted to a par- 
ticular problem of the quadrupolar ordering in CeBg. The theory naturally 
explains the growing specific heat jump at Tq{H) with increasing magnetic 
field. The role of the quadrupolar ordering in the formation of the magnetic 
ordering, as well as the possible critical experiments and applications to other 
rare-earth compounds, are discussed. 
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I. INTRODUCTION 



The aim of the present paper is to discuss the nature of the quadrupolar ordering in 
CeB 6 . This compound is classified as a dense Kondo system. With decreasing temperature, 
the resistivity grows logarithmically, attaining its maximum at T ps 3.2 K ffl. The Kondo 
temperature was initially estimated as T# ~ 8 K H. Later this value was significantly 
revised to a value of T# « 1 K from the experimental data on the magnetic susceptibility 
versus temperature This revision was caused by an unusual picture of the crystal field 
splitting, revealed in the Raman and neutron spectroscopic measurements ||. It is well- 
known that the crystalline electric field (CEF) of cubic symmetry (the elementary cell, 
containing the Ce ion with its boron environment, is shown in Fig. 1) results in splitting of 
the Ce 3+ multiplet (4/ 1 , J = 5/2, S = 1/2, L = 3) into a T 7 doublet and a T 8 quartet. The 
ground state of Ce 3+ in CeB 6 is realized as the well isolated T 8 quartet, and the T 8 — T 7 
CEF gap has been determined as 47 meV [|J . Prior to the results of Ref . [|J many difficulties 
in interpretation of experimental data had arisen in connection with incorrect assumptions 
on the multiplet splitting: The ground state level had generally been ascribed to T 7 (cf., 
however, HU). The quadrupolar and magnetic transitions were proved in the specific heat 
measurements 0-0, NMR [PUTT]] and neutron diffraction ||I~2"| , |i"3"f studies. Since the typical 



ordering temperatures (quadrupolar, Tq ps 3.3 K, and magnetic, T N ps 2.4 K) are much 
smaller than the CEF splitting, for low energy phenomena with T not exceeding several 
tens Kelvin, it should be legitimate to neglect a T 7 contribution and to deal with Tg only. 
The quadrupolar ordering is characterized by the following features: 

• There are two lines in the T — H phase diagram, which separate the anti-ferro- 
quadrupolar (AFQ) phase from the complex antiferro-magnetic (AFM) phases (see 
Fig. 2) and from the disordered (D) phase. The AFQ-D transition line, Tq(H), ex- 
hibits a highly anisotropic behavior. Starting at Tq « 3.3 K, Tq increases with H at 
not very high magnetic field, and increases linearly, dTQ(H)/dH > 0. The re-entrant 
behavior of Tq(H), predicted theoretically in [fLlp , has not yet been confirmed experi- 



mentally up to magnetic field of 18 Tesla [ 13), 15 ,|16|| (for recent experimental data, see 
Fig. 3). It is worth noting, that the estimates from below for the values of the critical 
field at the re-entrance and Tq = are 18 Tesla and 60 Tesla, respectively, according 



to Ref. ill 



There are contradictory AFQ patterns obtained in different microscopic measurements: 
neutron [0,[L5[], NMR |TT| and fiSR [17]]. The interpretation of neutron experiments 



is consistent with the AFQ patterns of the Q = [^hh] modulation, whereas the NMR 
and /xSR measurements display more complicated AFQ structures. Note, that until 
now the above-mentioned microscopic methods, i.e., 11 B-NMR, neutron diffraction, as 
well as /xSR, are used in non-zero magnetic fields which generate magnetically ordered 
states. Indeed, picking up an identical modulation with the AFQ state, a magnetic 
ordering is a secondary effect with respect to the primary quadrupolar ordering. The 
theoretical approach developed in Ref. selects the [m]-structure as energetically 
preferential. 

The specific heat jump at Tq appears to be of order of magnitude smaller than its 
counterpart at Tjv (H = 0) [|7||§- This points out an important role of fluctuations 



at the D-AFQ transition. This circumstance has been taken into account in [14 
by employing the spherical model description of the effective spin Hamiltonian. The 
specific heat jump on the D-AFQ transition line grows with H ||. 

The main features of the magnetic ordering have been presented in [Of : 



With a magnetic field applied along [111], the AFM-structure is characterized by the 
wave- vector, either ki = [|||] or k 2 = [|||] (the single-k structure at sufficiently high 
magnetic field), by a couple of k's, ki and k 2 (the doublle-k structure in moderate 
fields), and by a mixture of differently oriented domains at weak magnetic fields (see 
Fig. 2, where all these magnetic phases are sketched out). 

The Bragg peaks at the wave- vectors k x and/or k 2 are accompanied by ki = [ttO] 



and/or k2 = [77O], respectively [13]. Their occurrence is a sign of a crucial role of 



the AFQ modulation, Q = [220]' ^ n ^ ne formation of magnetic structures. A possible 



double-k structure, identified in [13], is shown in Fig. 4. 

In the next Section an effective "separation" of spin and orbital degrees of freedom 
is carried out by introducing a formalism, according to which magnetic and quadrupolar 
moments can be properly described by means of two Pauli matrices, cr and r. Section III 
concerns the analysis of the AFQ ground state and relevant excitations. In Section IV, on 
the basis of the two relevant interactions, Zeeman and quadrupolar, and using the spherical 
model for picking up these interactions, we are able to determine the shape of Tq(H). It 
occurs to be strongly anisotropic in the T—H plane. Despite a perfect cubic lattice symmetry, 
such a strong anisotropy is due to a spacial anisotropy of the quadrupolar interaction. Its 
conventional form, following from the Coulomb's interaction, gives rise to a very soft mode 
of r-excitations in the particular case of a simple cubic lattice. Experimentally, strong 
fluctuations are indicated by a small specific heat jump at the D-AFQ transition. This is a 
reason for employing the spherical model which is an appropriate tool for describing systems 
with developed fluctuations. 

The spherical model is applied for deriving analytical formulae for the specific heat 
near the AFQ-D transition. We also outline how the magnetically ordered state can be 
generated by the quadrupolar interaction via quantum fluctuations of orbital-like "spins", 
rs. In Section V the o — r representation is used for the case of a single / hole (configuration 
/ 13 ), which is likely ascribed to the rare-earth compound TmTe. In the concluding Section 
we discuss what kind of experiments could be critical for establishing the nature of the AFQ 
order unambigiously. 



II. THEORETICAL PREREQUISITE 
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A. Representation of moments through the Pauli matrices. Zeeman interaction 

We represent the set of the T 8 states with use of the \J Z ) (abbreviation for \L, S, J ', J z )) 
basis in the following form: 

^i,± = y||±5/2) + y||T3/2), ^,± = |±V2). (1) 

The quartet constituents in Eq.(|l]) are labelled in such a way in order to make use of the 
Pauli matrices, cr and r, convenient. For each £, the Kramers doublet ipe,± is defined as 

a z ipi± = ±l>4>£,± , <r+i/>t,- = ipe,+ (cr-^e,+ = i/>t,-) • (2) 

The orbital doublet ipi iCr and ip2,a can be suitably defined with using the pseudo-spin operator 
r as 

1 1 

7y01,<x = , T z 1p 2 ,a = ~2^ 2 ' a ' T+ ^ 2 > a = ^ ( T -^ 1 > CT = ^2, CT ) • (3) 

This representation was proposed in |J; however, expressions for the magnetic moment in 
terms of cr and r, given in ||, are oversimplified. 

In order to derive formulae for the moments (J, S, L, M), we need to calculate the 
matrix elements of, say, J over the set {ipe, a }'- 

<^l,±\Jz\^l,±>= ±^,<^2,±\Jz\^2,±> = ±L 

o 2 
2 2 

<^2, T \J±\^X,±>= -j=,<lkx\J±\ll>2,±>= -J=, (4) 

5 

<^1,±|J±|^1,=F> = ^<^2,±\J±\i'2, T >= 3. 

Within the Russell- Saunders scheme, the matrix elements of the moments can be obtained 
from their J counterpart in accordance with g-i actors of the f l multiplet: 

18 6 

<.. S.. >=--<.. J.. >, <..L..>= - <..J..>, <..M..>= -fM B <:j:> ■ 
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Using the matrix elements (|j) we can express the operator of magnetic moment M, which 
is associated with T 8 as follows: 

g 

Mi = 2/^B<7i(l + -Tj) , i = x,y,z, (5) 

where 

T - T - -- ^ T - ^ 

The derivation of formulae (||) from the set of matrix elements (f§) is outlined in Appendix 
A. Note, that the r^-component is not involved in (|5|). For the Zeeman interaction 

Hz = -H i J2M i (r), (7) 

r 

we shall use representation (||)-@. In (|7|) the sum runs over the Ce lattice sites. As usual, 
summation over repeated indices {i, Cartesian coordinates) is supposed. 

Let us discuss some simple properties of Hamiltonian (0), which are important in exper- 
imental applications to CeBg. 

1. If the "orbital", i.e., r-subsystem exhibits some AFQ order characterized by the mod- 
ulation vector Q, then at H ^ the effective Zeeman term, acting on as, produces, 
first, the uniform a-component, and second, the Q- modulated a- components. Both 
are absent in zero field. As a result, a uniform magnetic field causes the Q- modulated 
magnetization. This property of the AFQ phase has been used in neutron, NMR 
and /xSR experiments. In the weak-field region, with H not exceeding few Tesla, the 
Q-modulated magnetization is linear in H. 

2. In the magnetically ordered phase (see Fig. 2) the Bragg peaks are related either to the 
single k-structure (either k x = [|||], or k 2 = or to the double k-structure (ki 
and k 2 ). These peaks at k x and k 2 are accompanied by the Bragg peaks at k x = [ttO] 
and k 2 = [||0]. This fact can be easily understood if we note, that magnetization ([5]) 
is related not only to the a-modulations (wave vectors kx and/or k 2 ), but also to the 



(er ■ r)-modulations. The latter correspond to wave vectors ki + Q and/or k 2 + Q with 



3. Non-interacting r 8 ionic states can be realized practically, say, in La i _ :E Ce :E B 6 . Owing 
to a non-trivial form of the Zeeman interaction, magnetization is not aligned with H 
except for a few special if-orientations, e.g., [001], [110], [111], and their equivalents. 
At fixed H the bigger energy gain is achieved for directions of the [001]-type. This 
kind of H-anisotropy is an inherent property of the well-isolated T 8 states. 



Not only the vector moments, but the quadrupolar moment —x,y,z) on a Ce 

site as well, can be expressed in terms of the cr and r operators. For calculating the matrix 
elements of Qij over the set {ipe^}, 



we can employ the Wigner-Eckhart theorem, according to which these matrix elements are 
proportional to the operator equivalents: 



Given below are the matrix elements of the quadrupolar moment; we measure them in units 



Q = [HI1- 



B. Quadrupolar interaction 





of 



Qo =<^l,a\Qzz\^l,a>, 




(8) 



The Q matrix elements can be classified as a-independent 
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<lpl,a\Qzz\'lpl,<T>=-<?p2,*\Qzz\ , lp2,a>= 1, 
< tpl ,a | Qxx | ^1 ,<r > = < i>l ,a | Qyy | V>1 ,<t > = ~ \ , 
<^2,a\Qxx\i'2,a> = <^2, r 7\Qyy\lp2,a>=^, 
<?p2,a\Qxx\^l,a> = ~ < Qyy \lp2,a >= Y , 

as well as the cx-dependent 

<^2,+ \Qxy\^l,+ > = <^l,-\Qxy\^2,-> = i^, 

<lff2,-\Q V z\i>l,+ >=<^2,+ \Q Vg \^l,->=i^; ( 10 ) 
<^2,+ \Qzx\^l,->= ~ <lh t +\Q»x\'^i->= 3 §- 

We omit the Hermitian conjugated matrix elements in Eqs.(|9l)-(|10|). 

The matrix \\Q\\ can be written in the operator form as (see Appendix A for elementary 
explanations) : 



22"a; /i 2 fly 



\\Q\\ = Qo 



V 



(ii) 



l^y f^x ^ 

where /Xj = ^r y crj. The fact that contains the a- variables signals that the quadrupolar 
interaction can be responsible not only for pure orbital interactions, but also for magnetic 
interactions. 

The dependences of Mj and Qij on a and r determine the time-reversal properties of 
the cr and r components. It is evident from (|5]) that a — > — cr under t — > — t, whereas 
t x and r z are unchanged. The off-diagonal components of require r y — > —r y under the 
time-reversal transformation. 

We suppose that the predominant contribution to the interactions of Ce ions in CeBg 
comes from their quadrupolar interaction, the role of which in Ce compounds was first men- 
tioned by Bleaney ]Tj| (for a discussion on various forms of the quadrupolar interaction see, 
for instance, |I9| , |20|| ) . We accept the form of the quadrupolar interaction of the electrostatic 
origin, which is free of any model assumptions. Thus, our consideration is confined to the 
Zeeman and quadrupolar interactions: 
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^qd = H ^Aij,mn{r -r')Qii(r)Q m „(r / ) , (12) 

r^r' i...n 

where Aij tmn (r—r') is determined by the interaction V q (r — r') of two quadrupolar moments 
located at r and r'. The latter is given by 

V q (r) = —{2Q ij (0)Q ij (r) -20 g ij (0)Q iw (r)n i n m 

+35Qi j (0)Q m „(r)n i ra j n m ra n )}, m = Xi/r . (13) 

Thus, we obtain from (|13|) 



)} (14) 
Evident are the following properties of with respect to permutation of indices: 

J = 71 •■ A =71 

^ij,mn ^jt^n ^i^nm ^mrijjj- 

The diagonal elements of matrix ||<5|| give rise to the order parameter which is transformed 
according to representation T 3 characterized by two components (t x ,t z ) of r, while the 
off-diagonal elements are related to symmetry T 5 and are characterized by the vector ft. 
Keeping this in mind, we can rewrite 7i q d as follows: 

^qd = H [A a p{r-r')T a (r)Tp(r') + Bij(r-r')ni(r)nj(r')] , (15) 

where the Greek indices a, (3 prescribe summation over x- and z-components only. The 
expressions for A a $ and Bij are given in Appendix B. The Hamiltonian in Eq.(|15|) represents 
an evident separation of the orbital-like and spin-like parts. 

The magnetic exchange interactions are not relevant for a theoretical analysis of the 
AFQ ordering in CeB 6 . This applies to a major part of the phase diagram outside its 
low-temperature-and- weak- field part. The latter requires the RKKY- and Kondo-like inter- 
actions to be included. 
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III. TOWARDS QUADRUPOLAR ORDERING 



Taking the Fourier transform of the Hamiltonian (fl5|), we arrive at its k-diagonal form: 
n qA = E {^K, a r^ + B*MAO*kj} • (16) 

k 

We use the notation ^4 a ^[k] and f^fk] for the Fourier transformed coupling constants at 
general k. For high symmetry points of reciprocal space, such as [000], [±§§], [||0] and [00|], 
as well as along the cubic edge [||ac], the Fourier transformed Hamiltonian (Pj) becomes 
completely diagonal: 

^ = E {^kK.1 2 + At\r Kz \ 2 + 1/ikJ 2 ) + Skl^l 2 } • (17) 

k 

Table I shows the result of numerical calculations for the coefficients in units of Qq/ci 5 where 
a denotes the lattice constant. Let us estimate the order of such an energy unit. In doing 
so, we return to definition and then, performing the radial and angular integrations, we 
get: 

Qo = -^e <rf> . 



Then, the energy unit becomes 



16<r f >Ve» (ig) 



35 a ) a 



For CeB 6 , the lattice constant a ~ 4 A, the /-electron radius rf ~ 0.4 A, the lattice Coulomb 
unit e 2 /a ~ 3 eV, and we arrive at the Qo _un h of order 1 K (cf., however, (23). From Table 



I one can see that the coefficients are small as compared to the dominant ones, ,4s. 
Additional smallness of the i3-terms comes from the fact that the maximal value of [if is 
16/3 times smaller than the maximal value of t 2 . Thus, it seems appropriate to simplify the 
model by neglecting the £>-terms and to employ the simplified version of Ti^ in its purely 
orbital r-form: 

^orb = E E^/3( r - r/ )^(r)r /3 (r') . (19) 
10 



According to Table I, the global energy minimum could be achieved at k = Q. Not only 
the high-symmetry points of reciprocal space, but also wave vectors of a general position 
have been checked numerically in order to identify Q with the global energy minimum. It 
is necessary to emphasize that the energies at k = [||0] and Q are only slightly different. 
This is an indication of pronounced soft modes along the directions of cubic edges, i.e., 
[o k I] an d § < K < §■ Thus, competing AFQ patterns create fluctuations 

which should significantly decrease the AFQ-D transition temperature, as compared to the 
mean-field estimate. The wave vector Q is consistent with the AFQ patterns which have 
been found experimentally by the Grenoble group |13| . 

At the two points of reciprocal space, namely, and Q, we have A xx = A zz . Then the 
Fourier transform of 7Y or b takes a planar form 

^.k(T-k,xTk,x + T-k,2T k)Z ) . 

In other high-symmetry points, [||0] and [00^], 7i or b exhibits an easy-axis form with non- 
equal values of Af. and Af,. It is also valid for a general point of reciprocal space, but, in 
general, the off-diagonal component A xz is non-zero, and the easy-axis should be different 
from either x-axis or z-axis. 

It is worth noting, that searching for the ground state energy of the classical vector 
field r by using a Fourier transformation of Hamiltonian (|19D (also known as the Luttinger- 
Tissa method) would be a standard procedure, if the Hamiltonian were invariant under 
the homogeneous r rotations. Nevertheless, although the rotational symmetry of 7i rb at 
k = [||0] and k = [00~] is broken, the energy values listed in Table I are rigorous. 

A. Magnetic ordering due to electric quadrupolar interactions 

In this section we consider a quantum effect, namely the zero motion of the r "spins" 
with respect to the AFQ background. In fact, when a decoupling procedure is applied to 
the B terms in (|T5D, we obtain the effective spin-like Hamiltonian: 
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^m = ^E B&(r-0 <^(r)r y (r') > ^(r),^') = £ 4(r-i>,(r)a,(r') . (20) 

r^r' r^r' 

Tj, does not enter the Hamiltonian (|19D ; this is responsible for the formation of the orbital 
ordering. This would be a reason for neglecting all the contributions caused by r y , including 
interaction (|20|), were it not for quantum fluctuations of r. We put all the intermediate 
formulae, which determine our choice of the quantization axis, the spin-wave representation 
of r's, etc., into Appendix C. 

In the spin-wave approximation Hamiltonian (^) becomes 

n sw = J2 (^i(q)^q + i# a (q)(&- q 6 q + b\blS) , (2i) 



where 



K 2 (q) = ^(A zz [q]sm 2 4> + A xx [q]cos 2 ^-A xz [q\sm24>),q = q-Q, (22) 



and 



K 1 (q) = K 2 (q)-A Q . (23) 

For definition of 0, see Appendix C. The energy gain Eg), which occurs due to the zero- 
motion of r's is a straightforward result of the Hamiltonian (|2"ID diagonalization: 



2 

The correlation function 



WOlr^D-iE^Vf^f (25) 



appears to be non-zero although it decays exponentially with distance r. 

E^ is a periodic function of with periodicity n/3. In fact, using equations (|22|)-( 
and definitions of -4,'s given in Appendix B, one can rigorously prove that under the transfor- 
mation (q x , q y , q z ) — > (q z , q x , q y ) K\{q) and K 2 (q) remain unchanged, if is simultaneously 
shifted by n/3. 
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Numerical calculations show that <fi = 0, 7r/3, 27r/3, etc., are related to equivalent min- 
imums of E^. Using one of them, say, at <p = 0, we calculate the correlation functions 
( p5|) numerically. The sign of the first neighbor correlators is negative: (T y (0)T y (a x )) = 
(r y (0)T y (ay)) w -0.0293, (r y (0)r y (a 2 )) « -0.0058. This is a reminder of the AFQ ordering. 
Among the second neighbors only (r y (0)r 3/ (±a :E ± a y )) —0.0082 are not negligible, all 
others are much smaller. For the resulting coupling constants of (^) see Appendix B. 

This curious mechanism which, in principle, leads to the effective magnetic interactions 
(see (^0])), could be a reason for magnetic ordering at temperatures much smaller than 
Tq, because, first, the B coupling constants of Hamiltonian ( fL6|) are much less numerically 
than their A counterparts, and, second, the additional smallness comes from the quantum 
fluctuations of r y 's. The low-temperature magnetism in CeB 6 is unlikely to be described by 
such an unusual mechanism. Such a mechanism would come into play only when all other 
magnetic interactions (mainly via conductivity electrons) were very weak. 



IV. AF-QUADRUPOLAR DISORDER TRANSITION 

In this Section we consider CeB 6 near the AFQ-D transition. For this we employ, fol- 
lowing Ref. [14|, the spherical model which is applicable to a system with well-pronounced 



soft modes. The purpose of this section is to determine 

• the shape of the AFQ-D transition line in the T — H phase diagram; 

• the singularity of the specific heat along this line. 



A. Spherical model and AFQ-D transition line 

From the behavior of the specific heat anomaly (which is tiny in the weak magnetic 
field region) a strong short-range AFQ order should exist above Tq(H). A magnetic field 
suppresses the fluctuations and makes Tq(H) higher. In order to pick up these features, we 
go beyond the mean- field approximation for Hamiltonian Tiz + 7"4>rb- The first step in this 
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direction will be generalization of the spherical model for two spins, <r and r. For taking into 
account the quantum effects, we impose the constraints <cr 2 (r)>=3/4 and <r 2 (r)>=l/2. 
The latter would be equal to 3/4 were it not for redundancy of the r y variable. Note, that, 
as shown in Ref. ||22j| , the decoupling of fluctuations in the spin-1/2 Heisenberg model leads 
to the spherical model with the constraint <<r 2 (r) >=3/4. Now the partition function reads 

dtr(T) / dr{r) exp /?{A CT (3/4 - <r 2 (r)) + A T (l/2 - r 2 (r)) - (H orh + H z )} (26) 

, -oo J —CO 

where the spherical conditions 

dJ~ J d\ T = dTjd\ a = [T = -T In Z) , 

which control the constraints, should be satisfied by an appropriate choice of A T and A CT . 
Gaussian integration over cr(r) and r(r) in ( p6|) is straightforward, that allows us to derive 
the free energy T of the spherical model. 

For definiteness, we inspect the particular case of H||[001]: This orientation is expected 
to favor the re-entrance of the AFQ-D transition line at smaller H as compared to other 
orientations. The singularities on the AFQ-D transition line, as well as its shape, will be 
examined as temperature decreasing, i.e., from the side of the D-phase. 

Performing the routine calculations, which are given in Appendix D, we get (cf., 
Ml) 



4 2 \8J \ T +A -z 2 A,j 
T f d " k . ln & (27) 



where 



2 J (2vr) 3 (A t + A zz [k] - z)(X T + A xx [k\) - (A xz [k\Y 



8\ 2 (vbH) 2 



7J K 
Two equations, 

1 /7\ 2 z 2 T r d 3 k 2X T + A zz [k}+A xx [k}-z 



(28) 



2 \8J (X T +A -z) 2 2 J (2vr) 3 (A r + ^[k] - z)(X T + A xx [k]) - (A xz [k}y 
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(29) 



and 

3 /7\ 2 z /7\ 2 z 2 / z \ 3T 



2 + 



4 \8/ A<j V8/ X a (\ T +Ao- z) V A r +^ --2/ 2A a 

Z±f^ X T +A xx \k] 

~ 2 \J (27r) 3 (A r +^-[k]-z)(A r +^[k])-(^-[k]) 2 ' 1 J 

are valid for the D-phase and determine A CT and A T as functions of T and H, through which 

the physical quantities can be then expressed. The integrals in the r.h.s. of Eqs. (|29|) -(|30| ) 

are not of that simple form as the Watson integral, which appears in the spherical model 

treatment of the 3D ferromagnet, but let us call them generalized Watson integrals. The 

AFQ-D transition line corresponds to the values of A r at which the denominator of the 

generalized Watson integrals turns zero at k = Q. Taken at the critical value 

\[ c \z) = \A Q \ + z , (31) 

Eqs . (|29|) - (pOj) determine Tq(H). This line is depicted in Fig. 5. 

At zero magnetic field, A CT = 2Tq and Tq = (Wi(z = 0)) _1 , as seen from Eqs.(|30|) and 

(p9|), respectively. We denote the 1 st generalized Watson integral (see Eq.(p9|)) taken along 

the transition line by W\(z). It can be written now as 

w / \ _ f _fk 2\A Q \+A zz [k]+A^[k] + z 

1{ ) J (2n) 3 (\A Q \+A^[k]+z)(\A Q \+A zz [k])-(A^[k]y- 1 ) 

In order to estimate the importance of fluctuations, we can compare the spherical model 
and mean-field results for Tq. The mean-field approach yields Tq IF (H = 0) = \Aq\/2 = 5.37, 
and Tq(H) decreases monotonically with increasing H fL4j| . 

An important property of Tq(H) should be mentioned in connection with Fig. 5 (cf. 
Fig. 3), i.e., Tq(H) grows linearly with H at not very high magnetic fields. This feature 
is consistent with experimental findings. Mathematically, this behavior follows from the 
properties of the generalized Watson integrals: their expansions in a small parameter z 
are Cq — Cj\/z, where all four values of c are positive. Then, because in the leading 
order T Q (H) - T Q (0) « (Wi(0) - W x {z)) / '(^(O)) 2 , as it follows from Eq.@, the weak- 
field behavior of Tq becomes clearly understood. To complete this proof, let us represent 
(Wi(z) - Wi(0)) as follows: 
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W (0) W(z)-zf d ' k (lAjl+^W + ^W m 
WM - W 1 (z) - zj (27r)3 V ([k};z)V([k};z = 0) ' (33) 

which is positive. We denote the denominator of integral fl32|) by £>([k];z). At z — 0, 
integral in (|3^) becomes singular at small 5k = k — Q. In fact, the numerator behaves as 
5k 4 , whereas the denominator oc 5k 8 . The integral would diverge as |5k| -1 were it not for 
the cut-off at |5k| 2 ~ z. Thus we arrive at H\(0) — W\(z) oc sfz. 

It is worth noting, that the spherical model is a reasonable tool for picking up strong fluc- 
tuations. However if fluctuations are effectively suppressed, as is the case in high magnetic 
fields, the spherical model performs less satisfactorily. 



B. Specific heat 

Here we confine our consideration to the vicinity of the AFQ-D transition line and mag- 
netic fields weak as compared to Tq (hbH <C Tq). The equations which determine thermo- 
dynamic behavior are (|29|) -(|30l), where we can neglect all high order terms in H (H n with 
n > 2). Using the spherical conditions (<9jF /dX a = dT jdX T = 0), one can obtain for the 
entropy (S = —OT/dT): 

3 ttT 5 If d 3 k (ttT) 2 

" 2 n X a + 2 + 2J (2tt)3 11 (X T +A zz [k]-z)(X T +A™[k])-(A xz [k]y { } 

The specific heat at constant field, Ch = T (dS/dT) H , is now determined as follows 

5 Td\ T r d 3 k 2A T + A zz [k] + A xx [k] - z 



C 



H 



2 2 dT J (2tt)3 (A t + A zz [k] - z){\ T + A**[k\) - (A xz [k}) 2 
3TdX a Tz dX a r d 3 k X T + A xx [k] 



2 A (j dT 2A(j dT J (2tt) 3 (A t + A zz [k] - z)(X T + A**[k]) - {A xz [k}) 2 
With using Eqs.(E^)-(PD|) we can transform into the simple form: 



(35) 



5 lrfA T 3dX a 

CH -2~Yd 7 T~TdW (36) 

Let us consider a fixed magnetic field and T >Tq(H), i.e., T = Tq(H) + ST. At 5T <^Tq, 
X T = \Aq\ + z + SX T is supposed to be slightly different from Xf'(z) (see (PH)). In order to 
find out SX T (ST) we return to Eq. (|2~9|) which takes the following form at z <C Tq: 
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- = -(Tq(H) + 5T){W x {z) - ^8X T w(z/5X T ; 5X T ) ) , (37) 
where w is positive: 

w(z/oX T ; oX T ) 



V6X T 



d 3 k 


2|Aj 


+ *4 zz [k] + A xx [k] + 25A T + z 


(2vr) 3 (I^Q 


+A xx [k]+SX T + z){\A Q 


+A zz [k]+5X T )-A xz [k} 2 



(38) 

In Fig. 6 if is depicted as function of 5X T for three values of the ratio z/5X T . It can be 



remarkably well approximated to the form wo(5X T ) ■ (1 + a(5A T )y z/SXr)^ 1 , where w (x) = 
w(0;x) (see the upper curve in Fig. 6) and a(x) are both weakly dependent functions of 
the argument. For example, a(0) ~ 0.7 and a(l) ~ 1.15. Note, that a and wTq 2 are 
dimensionless quantities (w Tq 2 ~ 1). Thus, in the leading order, when 5X T <C z, which is 
equivalent to 5T <C HbH, we obtain: 



and, because Wi(z) ~ T Q 



-i 



^ a(0) Wi(z) /- 
w (0) T Q 



^ - a(0) 1 ^ (39) 



In accordance with Eq.(|3~0| 



dT w (0)T 2 



X <7 = 2T + 0(( f M B H) 2 /T Q ). (40) 



Therefore, in the leading order in H one can obtain from (B9), (0) and (|40[) : 



lim CB = 1 _^'i^. (41) 



Below Tq(H) X t does not depend anymore on 5T and remains equal to \Aq\ + 2. Thus, we 
get 

lim C H = 1 , (42) 

which, together with ([4lD , determines the specific heat jump, increasing linearly with i?. 
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At Tq S> 5T ^> fisH, i.e., beyond a narrow vicinity of the transition line, we arrive at 
Ch, decreasing linearly with ST: 



V. WELL-ISOLATED r 8 -LEVELS: A SINGLE F-HOLE CONFIGURATION 



Now we deal with the following quantum numbers: J = 7/2, 5 = 1/2 and L = 3. In the 
crystal field of cubic symmetry the eightfold multiplet splits into the T 8 quartet and two 
doublets, T 6 and T 7 . If the crystal field Hamiltonian allows the quartet to be realized as 
a well-isolated ground state level, then we confine our consideration to the following set of 
wave-functions: 



Listed below in units of hb are the non-zero matrix elements of M z , M + and M_: 



<V>1,+ 


M z \iJ 1A 


>= 


~<1pl 


-|M 2 |^i, 


>= 


44 
21 


<^2,+ 




>= 






>= 


4. 

7' 


<^2,- 


M + |Vi r 


->= 


<*I>1,+ 


M_|^ 2 ,_ 


>=- 


16 
7^ 


<^2,+ 


M_|Vi r 


>= 


<V , 1- 


M + |^2,+ 


>=- 


16 


<V'l,- 




r> = 


<V>1,+ 


M + |Vi,- 


>=- 


40. 

21' 


<^2,- 




h> = 


<V>2,+ 


M+1^2,- 


>=- 


24 
7 - 



All these matrix elements are in accordance with the operator expression: 




Zeeman and quadrupolar interactions 




(43) 



Mi = --/iWl + jTi) 



i = 



x,y,z. 



(45) 



Note, that the only difference between Eqs. (flop and (|5|) is the coefficient (-8/3 instead of 



2). 
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Now let us determine the matrix \\Q\\. For the unit, we take 

Qo = <^l,a\Qzz\4>l,*> ■ 

The diagonal components of the quadrupolar moment have the same non-zero matrix el- 
ements as in Eq.(|^), whereas the matrix elements of the off-diagonal components are six 
times larger than the corresponding elements of (|10|). Thus, the form of \\Q\\ is the same as 
in flriD , but the off-diagonal "vectors", which transform in accordance with T^, are defined 
now as /x = 3V3T y (T. This circumstance makes the problem of quadrupolar ordering less 
transparent as in the single /-electron case: The part of H q d (see Eq.([L^)), which reflects 
the fi — fi interactions, becomes as important as the r — r part. 

There are a few compounds of cubic symmetry, based on the single /-hole ions, which can 
be the candidates to realize the Tg quadrupolar ordering. Among them we would mention 
YbBi2 [^3| and TmTe p4 |. There is another important difference between these face-cubic 



centered compounds and CeB 6 (recall, that Ce ions are arranged in a simple cubic lattice). 
Namely, fee compounds do not exhibit such pronounced soft modes as in the simple cubic 
case. It is not our purpose to give a detailed analysis of the fee situation in the framework 
of the realistic quadrupolar interaction (|T^). We mention, however, that the ground state 
of the analog of H Q rb (see (|19"D) is of the AFQ-type, which is related to the high-symmetry 
points (Q x ,Qy,Qz) of the Brillouin zone boundaries (see Fig. 7): E(k = Qj) « —8.85 at 
them. The order parameter at, say k = Q^,, corresponds to the Ising-like symmetry with 
( T z) — and (t x ), altering the sign from layer to layer. The direction of low-lying excitations 
coincides with [001], but the mode is not a soft one: E(k = 0) — E(k = Q) « 3.99. 

VI. DISCUSSION AND CONCLUSIONS 

The problem of quadrupolar ordering in CeB 6 seems to be well-defined provided we 
confine our attention to the Zeeman energy and direct quadrupolar interactions. The un- 
usual form of these Zeeman and quadrupolar terms owes to the well- isolated r 8 quartet. 



19 



Instead of dealing with Stevens operators, it is more convenient to introduce the spin-like, 
cr, and orbital-like, r, operators. However, in the low-temperature and weak-field region 
CeBg undergoes the magnetic phase transition (see Fig. 2), which results in the appear- 
ance of complicated magnetic structures. For their explanation it is insufficient to restrict 
ourselves to the above-mentioned interactions only, but indirect interactions via conductiv- 
ity electrons start playing an important role. Magnetic domains of different orientations 
have been identified in neutron diffraction experiments [25| for magnetic fields applied along 



[111], [110] and [001]. However the interpretation of neutron experiments |]l3| , |i1)i|2~5| occurs 
to be contradictory to recent fiSK measurements [|l7j]. As for the AFQ ordered state, the 
neutron NMR results are still in disagreement: The triple-k structure has been proposed 



by Takigawa et al. [11], whereas in all the neutron experiments the Q modulation has been 



reported (Q^oo!!) - These two interpretations are mutually exclusive. Unfortunately, Ref. 



|flT| is a short paper with not many details, we mention a few of them to show a significant 
difference of the triple-k (qi = [|00], q2 = [0|0], q 3 = [00|]) and Q modulated structures. For 
the modulated component of magnetization m(r) the following equation has been proposed 



in P 



m(r) = (-lf mi + (-l) m m 2 + (-l) n m 3 , (46) 

where r=(£,m,n), and the polarization vectors depend on the magnetic field direction 
(H || (cos 0i, cos #2, cos 6*3)) through the equations: 

mi = m 1 (6 l i)(cos0i, — cos #2, — cos 6*3) 

m 2 = m 2 (0 2 )(- cos 0i, cos0 2 , - cos0 3 ) (47) 

m 3 = m 3(03)(~ COS 0i, — COS 2 , COS0 3 ) 

An important property of {mi} is that mj(7r/2) = 0. The form of Eqs.(fl6|)-(|4~7|) is not 
transparent, it is easier to illustrate them with a couple of examples. 

In the first example H||[001], which leads to 0i = 02 = 7r/2, hence, mi = rri2 = 0. This 
corresponds to a degeneracy of the triple-k structure, whose realization now is a single-k 
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structure. According to (f47|) the modulated component of magnetization is arranged as 
shown in Fig. 8a. Our theoretical description results in the arrangement shown in Fig. 8b. 

In the second example H||[110], and #3 = 7r/2, Q\ = 62 = 7r/4 lead to 772,3 = 0, mi = m^- 
This degeneracy corresponds to a double-k structure. Figs. 9a and 9b are the NMR and 
theoretical interpretations, respectively. Note, that the Fig. 9a pattern reproduces itself 
under any translation along [001], while ms of Fig. 9b alter the sign under the translation 
by the elementary lattice constant. 

The preliminary /iSR results jl7| contradict both neutron and NMR, measurements. 



A critical experiment seems to be not difficult to achieve. It is connected with the X- 
Ray structural measuments at zero magnetic field. The translational electric symmetry of 
the Ce-lattice is broken below Tq\ although all the electric charges of Ce ions are equal 
to each other, the modulated component of the quadrupolar moment should contribute to 
the Bragg reflections at the transferred wave vectors of the [||§]-type. This would be a 
weak effect, caused by only one electron of the total number Z. However, instead of an 
usual X-Ray technique it would be possible to use synchrotron facilities for finding the Q 
or triple-k modulated structure (or something different from these two). Note, that the 
zero field experiment is more instructive because it allows to avoid the secondary effects 
in non-zero fields, i.e. formations of magnetically modulated structures. To complete this 
X-Ray discussion we give the on-site form-factor operator which is calculated over the set 
{il>e t cr} and related to the / electron only: 

/(q) =<Jo> +2<U> p 4( z i) + =T*(r)(16<j 2 >-P2(cos0 q ) - 5<j' 4 >P 4 (^ q )) 

+ 773 TX ( r )( 8<j2>P 2 ( Z l) +< ^> P I (*q)) COS20 q , 2 q = COS# q . 

# q and q denote the spherical coordinates of q relative to the z-axis. P^iz) are asso- 
ciated Legendre polynomials. Detailed calculations in connection with a concrete X-Ray 
(synchrotron) experiment will be published elsewhere. 

The experimental technique which is associated with the so-called 3 rd order paramagnetic 
susceptibility |26| can be also used for probing the quadrupolar ordering in CeB 6 . However, 
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it cannot yield an information about the microscopic arrangement of quadrupolar moments. 
This method is based on the extraction of the i7 3 -terms from M(H), when H is small. For 
our problem, the average magnetization can be written as: 

M Q = --L1 E { (M a (ri ) M/3 (r 2 ) M M (r 3 ) M„ (r 4 ) ) 

ri...r4 

-^M a (v 1 )Mp{v 2 )){M^{v 1 )M v {v 2 ))}HpH^H v (48) 

For further transformations in (J4"8|), we imply that, first, all the magnetic fluctuations are 
much smaller as compared to the quadrupolar fluctuations in the AFQ phase (at least 
near Tq), second, the diagonal components of the quadrupolar moment are responsible for 
ordering below Tq. Then, using the Wigner-Eckhart theorem we can decouple Eq . ([48"D 
according to the following scheme: 

M a {v)M p {v') - (^)\ 2 B J a (r)J p (r)S r y -> QV|Q \T a (r) + J(J+l)/3)5 a p5 ry . 

and arrive at the following equation M a = + where 

= -|| E ^E((^( r )^(r')> " {T a {v))(T,{v')))H a Hl (49) 

m! 0) = -|| (£f^E E w.(a(r)^w^wxw) 

HX(r)^(r))(J,(r)X^ (50) 

Eq.(f49~D includes the irreducible correlators of r x and r z . The on-site irreducible correla- 
tor enters Eq. (|50|) . The average of four J7"'s can be reduced to a linear- in-r expectation 
value which disappears upon summation over r. The < J J > 2 terms of (|50|) result in the 
contribution oc (t) 2 , which should produce a kink in a dependence versus T at Tq. 

The AFQ-D transition line was recently measured in fields up to 18 Tesla [|TIJ. In spite 
of its tendency to re-enter, this curve still displays the monotonic Tq(H) behavior. An 
optimistic theoretical prediction is ~ 25-30 Tesla for the field at which the re-entrance 
could start and ~ 80 Tesla for the zero-temperature critical field. The current experimental 
facilities are enough to examine the field region around 25-30 Tesla. 
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Although many experimental results can be explained by the present theory (see also 



nj,27J), still there remain a few puzzling facts. Among them we would mention the AFQ-D 
transition line whose experimental shape is practically independent on the field orientation, 
[001], [110] or [111]. Probably, such behavior could be ascribed to an unusual anisotropy 
of the Zeeman energy. In fact, if we consider an isolated T 8 ion, its ground state energy 
depends on the magnetic field direction as follows: 



E gs = ~^y^65 + 4 V / 270« + n\ + n\) - 74, (51) 

that is -11, -9.81 and -9 (in units hbH/7) for orientations [001], [110] and [111], respectively. 
For magnetic field of a general orientation, the vector of magnetization in such a param- 
agnetic state does not follow the same direction. In this connection, the experiments with 
diluted compounds Lai_ a Ce x B6 could provide an important information, if the crystal field 
still favors the Tg ground state. 

From the theoretical point of view it should be interesting to understand the symmetry 
and a microscopic origin of the magnetic interactions which govern the properties of the 
system at low temperatures. 
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Appendix A 



In order to construct the operator expressions for M, as well as for Qij, we employ Table 
II. It shows the operator connection between all four possible states. Using matrix elements 
(||) in combination with Table II one can obtain for M x = (M + + M_)/2, for example: 

M x = ^|(4>/3(t_<7_ + t + ct + + r + a_ + r_a + ) + 10(1/2 + r z )(a + + a_) 

+18(1/2 - r z ){a + + <r_)) = 2^ B a x {\ + ^(VSr x - r z )) (52) 

Another example is for Q xy (see Eqs.(|10|)): 

Q xy = Qo (^(r_(l/2 + <r z ) + r+(l/2 - <r z ) 

-z ^(r+(l/2 + cr 2 ) + r_(l/2 - a,))") = ^a z r y = (53) 



Appendix B 

The original parameters ^jj, m n( r ) gi ye rise to ^4 Q/ g(r) and jBy(r) whose angular depen- 
dence is derived below in accordance with Eqs.( p]|L2l|L5|) : 

^5 5 7 

A =A +A +AA +2A -AA -AA =—(l-3n 2 ) 2 + -(l-3n 2 )-- 

•st-zz •^■xx,xx\^-yy,yy\^ L ^-zz,zz\^'^xx,yy ^<f^xx,zz ^^yy,zz 24 zl z) g' 

35 2 2 2 5 2 7 

A XX 3A XX<XX ~\~ ^Ayy^yy GA XX ^yy (^£ ^J/) (} ^^'2) g 7 

•^iz V^3 ( Axx, xx ~\~ *Ayy, yy 2A XX t zz 2Ayy^ zz ) V^3 (t^. ^y) ^ 24^ ^^^"^"g^ ' 

. 35 22^2^ /i/i 35 2 2 ^2 

= AAy z ,yz = ~~^ n y n z ~Q U x ~ g ' #f ~~ ^A ZXjZX — ~ n z ' f T' x + — H.^ — — , 

„ . 35 22^2^ k> /iyi /"t 2 i \ 

£*z.z ^Axy,xy ~Q^x^y g^ 2 g ' ^££/ AJ\y z ,xz g^as^y V < n>z J ' 

5 5 

^J/Z ^A ZX ^ X y -^TLyTL z {y TL x 1) , i3 2a; AA X y^y Z ~ T'ijjTlj; (7jly l) . 

To shorten the A and *B expressions we skip their r-dependence, i.e. factor r~ 5 . 

Listed below are a few coupling constants of the effective a — a Hamiltonian (|20|) : 
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B xx (a x ) = 0.0049, B yy (a x ) = B zz (a x ) = -0.0195 
B yy (a y ) = 0.0049, B xx {^) = B zz [a y ) = -0.0195 
B zz {a z ) = -0.0010, B xx (a z ) = B yy (a y ) = 0.0039 
B zz (a x ± a y ) = -0.0011, B xy (a x ± a y ) = ±0.0006 



Appendix C 

Because of the rotational symmetry of the Q-modulated state, we choose the quantization 
axis, n , as lying in the (r x , r z ) plane. One of the two perpendicular directions, n 1; is taken 
also in the (r x , r z ) plane. Let the third direction coincide with r y : 

n = (sin0, 0, cos</>) 
rix = (cos0, 0,— sin0) 
n 2 = ( 0, 1, ) 

In the spin-wave approximation 

r • n = e iQr (f - tfb 
T-ni = |e iQr (6 + 6t) 

r-n 2 = i(6-6t) 

Hence, to the same order we obtain 



r x = e 



iQr 



~ — n l sin y 



r, 



iQr 



(54) 



sin0 + |(6 + 6"!") cos ( 
f| — cos0 — |(6 + 6"!") sin< 

'»-s( 6 ~ 6t ) 

Now operators (|54|) can be used for extracting the "classical" and "spin-wave" parts of H or b : 
J2 (A zz (r - r')r z (r)r z (r') + A xx (r - r')r x (r)r x (r) + A xz (r - r')(r x .(r)r 2 (r') + r z {v)r x {v'))) 

r^r' 

= E(7 - h %) E e tQr (A zz (r) cos 2 + A.(r) sin 2 +A xz (r) sin 20) 

q r 

+7 E E e iQ V qr (A,(r) sin 2 + ^(r) cos 2 - A xz (r) sin20)(6_ q + 6+ )(6 q + &L q ) . (55) 

q r 
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The "classical", i.e., operator- independent, part of the r.h.s of Eg. ((55]) is invariant under 
the rotation of n in the (r 2 , r x ) plane (0-rotation). In fact, this part appears to be 0- 
independent: AM.q/4. 



Appendix D 

Let us first divide both fields, cr and r, into uniform and modulated parts: 

<T(r) = <7(°) + 5r(r), r(r) = r (0) +r(r). 

Then, the coefficients of the linear in a and r terms in the exponential of fl26|) must be put 
zero, and we obtain: 

a z (r) : -2\ a af) + 2fi B H(l + frf) = 0; Sr x (r) : erf = 0; ^(r) : <t<°> = 0; 

f,(r) : -2(A T + ^ )ri°) + f /i^af = 0; f x (r) : = 0; f„(r) : rf = 0. 
This yields 

(0) _ 7 Z (o) /^S-H" + Aq 



(7: 



8 A T + v4o — ^ Ao- A t + Ao — z 

where z = (8fx B H/7) 2 / X a . 

Then, extracting the constants and the contribution of the uniform components of o- 

and r-fields (Zq), we arrive at the Gaussian integration over a and r in Z\ {Z = ZqZ{): 
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/oo roc ( \fo 

da{r) dT{r)exp(3 -A ff cr 2 (r) + —^ B Ha z {v)f z {v) 
, -oo J — OO L f 

-A r f 2 (r) - " r')f a (r)^(r')| 

where Z = exp(— NFq/T) and 



(56) 



-* = \^+®\> + T^)- (57) 

The a integration can be easily performed. It results in the contribution 

-^ = §Tln^, (58) 

and transforms A zz (r — r') — > A zz (r — r') — 5 r)I ./2;. 
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Thus, we can rewrite Z\ as exp(— N{Tx + Fi'^/T), and 
NJ 7 ^ r°° 

exp — - - Yl I rfT(r)exp-/3^^ Q/3 (r-r , )r Q (r)r / 3(r / ), (59) 

1 r J-oo r , 

where 

A zz {r) = 5 r , (A r - z) + A«(r), A»(r) = S rfi X T + A xx (r), A xz (r) = A xz (r). 
The functional integration in ( |55| ) can be straightforwardly performed: 

2 iV k ^[k]^[k] - (^[k]) 2 V 7 



27 



REFERENCES 

* On leave from Landau Institute for Theoretical Physics, Chernogolovka, Moscow Dis- 
trict 142432, Russia. 

[1] K. Winzer, and W. Flesch, J. Physique, 39, C6-832 (1978). 

[2] A. Takase, K. Kojima, T. Komatsubara, and T. Kasuya, Solid State Commun., 36, 461 
(1980). 

[3] N. Sato, S. Kunii, I. Oguro, T. Komatsubara, and T. Kasuya, J. Phys. Soc. Japan 53, 
3967 (1984). 

[4] E. Zirngiebl, B. Hillebrands, S. Blumenroder, G. Giintherodt, M. Loewenhaupt, J. Ca- 
rpenter, K. Winzer, and Z. Fisk, Phys. Rev. 30, 4052 (1984). 

[5] P. Burlet, J.X. Boucherle, J. Rossat-Mignod, J.W. Cable, W.C. Koehler, S. Kunii, and 
T. Kasuya, J. Physique 43, C7-273 (1982). 

[6] F.J. Ohkawa, J. Phys. Soc. Jpn. 52, 3897 (1983). 

[7] K.N. Lee, and B. Bell, Phys. Rev. B 6, 1032 (1972). 

[8] T. Fujita, M. Suzuki, T. Komatsubara, S. Kunii, T. Kasuya and T. Ohtsuka, Solid State 
Commun. 35, 569 (1980) 

[9] Y. Peysson, C. Ayache, J. Rossat-Mignod, S. Kunii, and T. Kasuya, J. Physique 47, 
113 (1986). 

[10] M. Kawakami, K. Mizuno, S. Kunii, T. Kasuya, H. Enokiya, and K. Kume, J. Magn. 
Magn. Mater. 30, 201 (1982). 

[11] M. Takigawa, H. Yasuoka, T. Tanaka, and Y. Ishizawa, J. Phys. Soc. Jpn. 52, 728 
(1983). 

[12] S. Horn, F. Steglich, M. Loewenhaupt, H. Scheuer, W. Flesch, and K. Winzer, Z. Phys. 

28 



B 42, 125 (1981). 

[13] J.M. Effantin, J. Rossat-Mignod, P. Burlet, H. Bartholin, S. Kunii, and T. Kasuya, J. 
Magn. Magn. Mater. 47-48, 145 (1985). 



[14 
[15 

[16 

[17; 

[18 
[19 

[20 
[21 



[22 
[23 

[24 
[25 



G. Uimin, Y. Kuramoto, and N. Fukushima, Solid State Commun. 97, 595 (1996). 

W.A.C. Erkelens, L.P. Regnault, P. Burlet, J. Rossat-Mignod, S. Kunii, and T. Kasuya, 
J. Magn. Magn. Mater. 63-64, 61 (1987). 

A. Lacerda and M. Torikachvili, unpublished. 

R. Feyerherm, A. Amato, F.N. Gygax, A. Schenck, Y. Onuki, and N. Sato, Physica B 
194-196, 357 (1994); J. Magn. Magn. Mater. 140-144 1175 (1995). 

B. Bleaney, Proc. Phys. Soc. (London) 77, 113 (1961). 

R.J. Birgeneau, M.T. Hutchings, J.M. Baker, and J.D. Riley, J. Appl. Phys. 40, 1070 
(1969). 

P.M. Levy, P. Morin, and D. Schmitt, Phys. Rev. Lett. 42, 1417 (1979). 



According to [f20] the value of the direct quadrupolar interaction in some intermetallic 
rare-earth compounds is of order 10~ 3 K. That estimate should be consistent with a 
large screening effect on a quadrupolar moment by outer electrons, which is unlikely the 
case of CeB 6 . 

J. Hubbard, J.Phys.C 4, 53 (1971). 

K. Sugiyama, F. Iga, M. Kasaya, T. Kasuya, and M. Date, J. Phys. Soc. Jpn. 57, 3946 
(1988). 

T. Matsumura et ai, unpublished. 

For a detailed description of orientations and periodicities of such domains we may 
recommend the following thesis: J.M. Effantin, Universite de Grenoble (1985); W.A.C. 



29 



Erkelens, University of Leiden (1987); A. Bouvet, Universite J. Fourier, Grenoble (1993). 
[26] P. Morin, and D. Schmitt, Phys. Lett. A 73, 67 (1979). 
[27] G. Uimin, Phys. Lett. A 215, 97 (1996). 



30 



FIGURES 





Temperature (K) 



FIG. 2 The low field part of the phase diagram. Positions of the lines, confining the magnetically 
ordered phases, depend on the orientation of magnetic field. 
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FIG. 3 The boundary between two phases, AFQ and D, as determind experimentally [|l6| by 
transport (magnetoresistance) and magnetic measurements, full and open circles, respectively. A 
line is included as a guide to the eye. 
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FIG. 4 One of the possible arrangements of magnetic moments in the double-k structure. 
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FIG. 5 The line of the AFQ-D phase transition according to Eqs.(26)-(27). Units for T and H 
are discussed in the text. 
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FIG. 7 The first Brillouin zone of reciprocal space of the fee structure. 
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FIG. 8 Arrangement of staggered magnetization for magnetic field applied along [001]; (a) NMR 
interpretation [11], (b) Q-modulated structure. 
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FIG. 9 Arrangement of staggered magnetization for magnetic field applied along [110]; (a) NMR 
interpretation [11], (b) Q-modulated structure. 
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TABLES 



k 


K 


Al 


Bt 


Bi 


111 

2 2 2 


-10.736 


-10.736 


1.789 


1.789 


2 2 U 


-10.478 


3.484 


-0.581 


2.909 


oo± 


2.139 


7.349 


-0.357 


-1.659 


000 


9.325 


9.325 


-1.554 


-1.554 



TABLE I. Strength of the quadrupolar interaction at high symmetry points 





IL+) 


IL-) 


|2,+) 


|2,-> 


IL+) 


{k + T z ){\ + a z ) 


(3 + r 2 )a_ 




r_(j_ 


IL-) 


(5 + r z )a + 


(£+r*)(±-<7*) 






|2,+> 




T + cr_ 


(\-T z ){\ + a z ) 


(5 - T 2 )cr_ 


|2,-> 


T+CT+ 


r+(5 - <r z ) 


(3 - r z )<7+ 


(| - t z ){\ - a z ) 



TABLE II. The operator forms for transformations \i,a) — > |^',cr'). 
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